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Abstract 

Several preliminary QCD results from e+e~ interactions at LEP are 
reported. These include studies of event shape variables, which are used to 
determine Qs and for studies of the validity of power corrections. Further, 
a study of color reconnection effects in 3-jet Z decays is reported. 



1 CKs from event shape variables 

Variables which quantify in some way the distribution of the particles of an event 
in momentum space, known as shape variables, are sensitive to the amount of 
hard gluon emission in the event. The distributions of such variables are there- 
fore sensitive to the value of the strong coupling constant, ttg- The distributions 
expected by QCD can be fit to the data distributions in order to measure as- 

To be useful for this purpose, the variables should be infrared and coUinear 
safe and insensitive to the electroweak physics which produces the event. Exam- 
ples of such variables are r = 1 — T, where T is the thrust; the scaled heavy jet 
mass, Mh; the total and wide jet broadening, _Bt and -Bw; the C-parameter, C; 
and the jet resolution parameter, ?/23, which is the value of j/cut in the Durham 
algorithm at which the event classification changes from 2-jet to 3-jet. 

All four LEP collaborations have measured these distributions at various 
center-of-mass energies and used them to determine as- The LEP QCD Working 
Group has performed a preliminary simultaneous fit to all of these distributions 
[1], which is described here. Such a combined fit allows a consistent treatment of 
the theoretical predictions and uncertainties, as well as of correlations between 
variables and energies. 

Each experiment measures the shape variable distributions and corrects 
them for detector resolution and acceptance, background, initial state radia- 
tion, etc. The theory predictions are calculated and, since these predictions 
are at parton level, corrected for hadronization using a parton shower Monte 
Carlo (MC) program such as pythia, herwig, Ariadne. The corrected theory 
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predictions are then fit to the corrected experimental distributions to determine 
To C(Q!g) the distribution of shape variable y is given by 



/pert(y;a;s) = 



CTtot dy 2ir 



1 da as 




where /3o = ^^5^) with rif the number of active flavors, and /z = x^-^/i is 
the rcnormalization scale, providing a parameter to use to vary the scale. 
Integration of the ERT matrix elements gives the values of A and B. This 
describes the data well in the multi-jet region, but not in the 2-jet region, which 
corresponds to small values of where emission of softer gluons is important. 
They can be included by summing, to all orders of Ug, the leading and next- 
to-leading order logarithmic terms in the expansion of R{y; as) = /pert dy in 
terms of L = ln(l/y). The two calculations can be combined if one is careful 
to avoid double counting. This is advantageous since it allows use of a fit 
range extending into the 2-jet region. However, it is not without theoretical 
uncertainty. There are two "matching schemes" to do this, the Log-R and the 
Modified Log-R schemes. The latter forces R to vanish above the kinematical 
maximum value of y by replacing L by 



where the parameters p and xl allow variation of the incorporation of the kine- 
matical limit. 

The data samples available cover three center-of-mass energy ranges: 20- 
85GeV from radiative Z decays where the radiated photon is removed from 
the event; Mz; and 133-206 GeV. At present only L3 measurements are used 
from the lowest energy range, although OPAL has recently released preliminary 
measurements [2]. The values of ttg from this new OPAL analysis agree well 
with the current world average and are therefore not expected to have a large 
effect on the combination. 

At present 6 shape variables at 14 energies are used giving 194 measurements 
in total, not all variables being available for all energies/experiments. To per- 
form the fit, the covariance matrix (194 x 194) of these measurements is needed. 
It is composed of four contributions: Vij = V^f** -|- ^ + V^^'"^ + V^>. The 
first term, the statistical uncertainty, is the easiest to evaluate. It is certainly 
uncorrelated between experiments and between energies. The second term is 
the systematic uncertainty in the experimental measurement. It is uncorrelated 
between experiments. Correlations within an experiment are taken as the "min- 
imal overlap" , V^^^^ = mm{V^^^, Vj^^). The third term is the systematic uncer- 
tainty in the hadronization correction of the perturbative calculations. While 
one might expect this to be correlated between experiments, since they all use 
the same programs, it seems that these correlations are small. This is presumed 
to be due to the fact that each experiment uses a different tuning of MC param- 
eters. These uncertainties are estimated from a comparison of the corrections 
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found using PYTHIA, HERWIG, ARIADNE. Large fluctuations arc seen in the val- 
ues from energy to energy, presumably arising from statistical limitations in the 
Monte Carlo samples. Accordingly, they were smoothed assuming a. 1/Q de- 
pendence, as suggested by power corrections, and all correlations were assiimed 
to be zero. The fourth term is the uncertainty in the perturbative prediction 
due to neglect of higher orders. There are certainly very large correlations not 
only between experiments, but also between energies and shape variables. How- 
ever, attempts to fit including such large correlations lead to highly unstable 
results. Accordingly, all correlations were set to zero. To evaluate the diagonal 
elements, the following procedure was followed. First, a nominal value of as 
was chosen, and the shape variable distributions were calculated for the default 
values = xi, = p = 1 using the modified Log-R scheme, giving the nominal 
distribution. Distributions were then calculated varying the parameters x^, .tl 
and p as well as using the Log-R scheme. The difference of these distributions 
with the nominal one is shown for a typical case [3] in Fig. 1. Then the value 
of as is varied using x^ = .xl = p = 1 to determine values of as which produce 
changes in the distribution as large as the largest differences found within the 
fit range. The difference between this value of ag and the nominal value is taken 
as the uncertainty on ctg. 

Using the above-described covariance matrix a least squares fit is performed. 
Because of the assumptions made on lack of correlations in the hadronization 
and theory uncertainties, the uncertainties found by the fit can not be expected 
to be correct. To determine the final hadronizaton uncertainty, the combination 
is performed 3 times, once for each of the MC programs. The result for ag is 
found using PYTHIA, the uncertainty from the rms of the results. The theory 
uncertainty is found by repeating the fit twice, using for all points asi ± V^it- 
The theory uncertainty is taken as half the difference. The result is 

a,(Mz) = 0.1201 ± 0.0003 (stat) ± 0.0009 (exp) ± 0.0009 (had) ± 0.0047 (th) 

Note that the theory uncertainty dominates. The results from different energies 
and shape variables are consistent, as shown in Fig. 2. 

2 Power corrections 

The power correction ansatz parametrizes the unknown behavior of as below an 
infrared matching scale, /ii, by an average, ao = — Jq^ as{k)dk. This leads to 
a power term, P oc l/^/s which shifts distributions of shape variables: f{y) = 
/pert(j/ — CyP) and increases their moments: (y) — {y)pen + CyP and (j/^) = 
(2/^)pert + 2(y)pert CyP. The factor Cy is a known factor, different for each shape 
variable, but P is supposed to be universal. 

DELPHI [4, 5] has analyzed both the distributions of a number of shape 
variables and their first moments. L3 [6] has analyzed both first and second 
moments. The first moments of different shape variables result in consistent 
values of ag, but differences of around 20% are observed in the values of ao- 
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Figure 1: Variations to determine the theoretical uncertainty on as- 
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Figure 2: Results for as at different energies and using different shape variables. 
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The situation is much worse in the analysis of the distributions, as shown in 
Fig. 3, and of the second moments. 

While one would have hoped that power corrections could have been used in- 
stead of Monte Carlo models for the hadronization corrections of shape variables, 
the results indicate that power corrections provide only a semi-quantitative de- 
scription. 
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Figure 3: Results of analysis of shape-variable distributions using a power cor- 
rection ansatz. 



6 



3 Color reconnect ion in 3-jet Z decays 

So-called rapidity gap events have been observed in ep and pp events and are 
attributed to color-singlet exchange. OPAL has investigated [7] whether such 
effects exist in e^e^ events by studying directly the distribution of particles 
within the gluon jet of 3-jet Z decays. L3 has taken a different approach [8]. 

The occurence of a eolor-singlet exchange within the gluon jet breaks the 
jet into two pieces, one which itself is a color singlet, and another which forms 
a color singlet together with the quark and antiquark. This has an effect on 
the color flow between jets. To quantify this, asymmetries are constructed 
comparing the flow between q and g with that between q and q, eg, = 
{—Bi2 + B23 + B31)/ (i?i2 + B23 + -B31), where Bij is the smallest angle between 
two adjacent particles in the region between jets i and 7, excluding particles 
within cones of 15° about the jet axes. A "Mercedes" topology is required with 
jet 3 identified as the gluon jet by a b-tag for jets 1 and 2 and an anti-b-tag for 
jet 3, resulting in a sample of 2668 events with a gluon jet purity of about 78% 

The asymmetry distribution is shown in Fig. 4 and compared to the expec- 
tations of various Monte Carlo models, both without and with a "color recon- 
nection" (CR) algorithm. Both jetset and ariadne without CR agree well 
with the data. However, when the GAL model of Rathsman [9] is included in 
JETSET or when the CR model in ARIADNE is used, the models disagree with 
the data. The agreement between herwig and the data is very poor both with 
and without its CR model. The failure of the CR models here suggests that 
they are also inapplicable to the case of CR in e+e~ W+W~ — > qqqq. 
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Figure 4: A color-flow asymmetry variable for 3-jet events compared to MC 
predictions with and without color reconnection models. 
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